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We generalize our previous linear result [T] in obtaining gravitational waves from 
our piecewise flat model for gravity in 3+1 dimensions to exact piecewise flat config- 
urations describing exact planar gravitational waves. We show explicitly how to con- 
struct a piecewise flat spacetime that describes an impulsive plane wavefront. From 
these wavefronts more general plane waves may be constructed. Further remarks 
are made on how this construction may be extended to non-plane (e.g. cylindrical) 
waves. 



I. INTRODUCTION 

In previous articles (2j [3] 't Hooft and the author have introduced a piecewise flat model 
for gravity in 3+1 dimensions, motivated by general relativity in 2+1 dimensions. In 2+1 
dimensions, the Riemann tensor is equivalent to the Einstein tensor. Consequently, the 
Einstein equation completely and locally determines the geometry in terms of the matter 
content. In particular, empty space is flat and has no local purely gravitational degrees of 
freedom, facilitating quantization of pure gravity [U [5] . A system of point particles generates 
a piecewise flat geometry, with the particles represented by moving conical curvature defects 
[5] . Such a configuration has locally only a finite number of degrees of freedom and may be 
quantized [H |S] . 

One obstruction to quantizing gravity in 3+1 dimensions is the presence of local gravita- 
tional degrees of freedom in the form of a spin-2 graviton field, which is non-renormalizable 
when quantized perturbatively. The idea behind the model discussed in [21 E] is to take a 
page from gravity in 2+1 dimensions, and eliminate all local degrees of freedom by enforcing 
that empty space should be completely flat. Local degrees of freedom are introduced as con- 
ical curvature defects. If we insist that the Einstein equation be satisfied, then the flatness of 
empty space implies that these defects are flat co-dimension 2 surfaces. Since these defects 
are interpreted as physical degrees of freedom, causality requires that only timelike surfaces 
appear. This distinguishes the model from other piecewise flat approaches to gravity such 
as Regge calculus [HI \10\ ■ 

The additional assumption that empty space is flat appears quite limiting. In partic- 
ular, at first sight it appears that it completely eliminates the possibility of gravitational 
waves. However, in the model the propagating line defects represent both the matter and 
gravitational degrees of freedom. For this reason the model allows both defects that con- 
tribute positively and negatively to the curvature]]] The appearance of positive and negative 
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1 It the defects are interpreted purely as matter then the negative curvature defects would correspond to 
matter with negative energy density, and should not appear. 
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curvature defects opens the possibility of geometries for which the Einstein curvature van- 
ishes when averaged over larger scales, while the Riemann curvature does not. Some such 
geometries could be interpreted as emergent gravitational waves. 

In our previous article pQ , we studied this possibility in the linearized limit of our piecewise 
flat gravity model. The two main results were that: 1) Configurations of moving defects can 
produce the energy-momentum distribution of an arbitrary dust cloud of point particles, 
despite the spacial extent of the defect lines; 2) There exists a large family of linear defect 
configurations with vanishing energy-momentum, whose metric perturbation corresponds 
with linear gravitational waves in general relativity. At the linear level, the model is able to 
reproduce gravitational waves as an emergent feature at larger scales. 

In this article we extend this result to exact gravitational waves. In contrast to linearized 
general relativity, it is not possible to erive completely general exact vacuum solutions of 
general relativity. Deriving exact solutions requires some form of additional symmetry. For 
example, it is known how to construct completely general gravitational waves with cylindrical 
symmetry 1 1 1 j . and how to construct exact gravitational plane waves [T2] . 

Here we will focus on plane gravitational waves and show how to construct a family of 
increasingly fine discrete configurations of defects that approach an exact gravitational wave 
in the limit that the configuration becomes continuous. 

Our strategy will be to first consider gravitational plane waves in the linearized limit, and 
examine what linearized distributions of defects approximate these waves (section [n]). Based 
on the linearized configuration, we will make an ansatz for an exact piecewise flat planar 
wavefront, which reduces to the linearized distribution in the linear limit. We show that 
this ansatz is indeed a bona-fide piecewise flat configuration, and that in the limit where the 
configuration becomes continuous it produces an impulsive gravitational wave solution of 
general relativity (section III). In section IV we combine trains of these wavefronts to form 
general exact gravitational plane wave solutions. Finally, in section [V] we comment on the 
construction of non-planar waves and the additional subtleties that those entail. 



II. LINEAR PLANE WAVES 

For simplicity we shall consider, without loss of generality, plane waves travelling in the 
positive ^-direction. Expressed in lightcone coordinates (u, v , x, y)Q such a wave takes the 
following form in linearized general relativity, 
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(1) 



where h^ u is the perturbation of the metric on the Minkowski background. In our previous 
article p] we saw how to produce such a metric perturbation in linearized piecewise flat 
gravity. The basic building block was a configuration we called a "laminar plane wave of 
defects" . A laminar plane wave of defects is a space-filling configuration of parallel conical 
curvature defects (i.e. line defects that have the same orientation and velocity), whose 
energy density varies as the amplitude of a plane wave with a wave vector perpendicular to 

2 In this article, spacetime metrics have signature ( — h ++), lightcone coordinates are defined as u = 
and v = ^m, and we use units such that c = Ti = 8ttG = 1. 
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Figure 1. A laminar plane wave consists of parallel line defects travelling in the same direction, 
whose densities vary as a plane wave in the direction of propagation. 



the defects. That is, given a point in our spacetime there is exactly one defect (spanned by 
4-vectors eP and w M ) in the laminar plane wave configuration that passes through this point, 
because all the defects in the laminar plane wave are parallel. The linear energy density of 
this defect is given by Oe 1 ^'" with some fixed wave vector k^, which is perpendicular to d^ 
and m m . 

In the linearized limit of our piecewise flat model an arbitrary configuration of defects 
may be written as a linear combination of laminar plane waves. In [1] we found that a 
gravitational wave of the form ([I]) was a produced as a superposition of laminar plane waves 
of defects with the velocity = (1, 0, 0, 1) and wave vector /c M = (u, 0, 0, u>). The orientation 
dp of the defects in these laminar plane waves can then be parametrized by a single angle, 
0, which we choose to be the angle between the defect lines and the x-axis. 

To be a gravitational wave the total energy density amplitude of the superposition of 
laminar plane waves must vanish. That is, if the energy density amplitude of the laminar 
plane wave with orientation is called 0(0), then this condition can be expressed as 



O(0)d0 = 0. 



(2) 



It was established in [T] that this is the only condition on the superposition of laminar plane 
waves of defects required for it to describe a gravitational plane wave. The polarization of 
the resulting gravitational wave was shown to be given by 



h + = / 0(0) cos(20)d0, 
Jo 

/*7T 

h x = 0(0) sin(20)d0. 
Jo 



(3) 



All higher moments of 0(0) are unconstrained. That is, there is a large family of linear 
configurations of defects that produce the same linear gravitational wave. We therefore have 
quite a bit of playing room for finding a linearized configuration that can be interpreted as 
the limit of some exact configuration of defects. 
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Figure 2. A wavefront of perpendicular positive (white) and negative (black) energy line defects. 



Geometrically, it would be convenient if fi(0) has as few discrete components as possible. 
Condition (|2]) implies that we need at least two components. This is, in fact, sufficient to 
create waves of any polarization. If we take 

ft(0) = a (5(0 - O ) - 5(<P - 0o - f )) , (4) 

then the polarizations are given as h + = 2acos20o and h x = 2asin20o. We thus have 
perpendicular components of opposite energy, with the amplitude of the defects determining 
the amplitude of the wave, and the orientation determining the polarization. 



III. PLANE WAVEFRONTS 

This linear configuration prompts us to make an ansatz for an exact configuration of 
defects. Further discretizing the linear distribution, we guess that a configuration consisting 
of subsequent grids of perpendicular positive and negative energy massless defects should 
suffice to approximate a gravitational wave. To analyse this configuration we first consider 
a wavefront located at u = ^£ = consisting of a single grid (see figure 2 ) . 

We first need to check if this ansatz defines a proper configuration of defects. In [2] 
the conditions that should be met at a junction of multiple defects were discussed. These 
can be phrased in terms of the holonomies of loops around the defect lines. They can be 
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summarized as saying that the product of the holonomies of a set of loops whose concatena- 
tion is contractible should be the identity. For the proposed grid configuration the junction 
conditions must be satisfied at each vertex of the grid. 

Since all the horizontal line defects in figure [2]have the same density, the holonomy of each 
simple loop around a section of a horizontal line defect should be the same, and similarly 
for the holonomies of simple loops around vertical line defects. If we call the holonomy of a 
simple loop around a section of a horizontal grid lines Q a and the holonomies of the vertical 
grid lines Qb, then at each vertex the junction condition is 

QaQbQa^Qb^ = I- (5) 

The holonomy of a lightlike defect is a null rotation, i.e. a Lorentz transformation that 
leaves a lightlike and a perpendicular spacelike direction invariant. In the PSL(2, C) rep- 
resentation of the Lorentz group, the holonomy of a general lightlike defect moving in the 
^-direction is represented as 

Qz = (I f) , (6) 

where ( is a complex number. The argument of ( gives the direction of the defect, while 
its modulus gives the energy density of the defect, i.e. the parabolic angle of the null 
rotation. The junction condition for n lightlike defects moving in the z direction with 
complex parameters 0, . . . , ( n therefore reduces to 

Cl + • • • + Cn = 0. (7) 

Furthermore if Q a has complex parameter ( a then Q a ~ l has complex parameter —( a . Hence 
the junction condition for the vertices of the grid ^ reduces to 

(a + Cb ~ (a ~ Cb = 0, (8) 

which is true for any ( a and Cb- Our ansatz configuration of a grid of perpendicular lightlike 
defects is therefore a bona-fide configuration of defects for any value of the amplitudes of its 
components. In particular, it is valid for the particular case we are interested in, where the 
components have opposite amplitudes. 

Our next step is to compare this configuration with an exact plane wave solution of the 
vacuum Einstein equation. For this we need to find a metric to describe the geometry of 
the grid of defects. 

Away from the grid the spacetime is flat. Hence, if we position the grid at u — 0, we can 
choose the metric for u < to have the familiar Minkowski form, 

ds 2 = 2&u&v + da; 2 + dy 2 . (9) 

For simplicity we choose to align the grid with the x and y-axesj^] such that the positive 
energy defects are aligned in the x-direction and sit at constant values of y = e{n + ^) for 
integer values of n and grid spacing e, and the negative energy defects are oriented in the 
y-direction and sit at constant values of x = e{m + ^). 

The defect and surplus angles of the grid lines can be oriented along the direction of 
propagation of the grid, i.e. the positive u direction. The cuts of the defect /surplus angles 



3 This should correspond to wavefront with "+" polarization, according to the linear analysis. 
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Figure 3. The geometry of one "tube" of spacetime.The deficit angles of the positive energy defects 
(white), cause the geometry to be squeezed in the y-direction as u increases. Similarly the surplus 
angles of the negative energy defects (black) cause the geometry to be elongated in the x-direction. 



then divide the u > side of the grid in square "tubes" of spacetime (see figure [3]). The 
Minkowski metric from u < side of the grid can be continued into each tube. The constant 
u and v slices in each tube start out as a square at the base (u = 0) of the tube. As u 
increases the square is elongated in the x direction and compressed in the y direction. 

We now construct coordinates on each tube, which can be patched together to one coor- 
dinate system on the u > side of the spacetime. The space 



T 



R> xRx [e(m - §), e(m + §)] x [e(n - §), e(n + §)] 



can be mapped to the tube behind the face of the grid centered on x = em and y 
the mapping, 

u i->- u, 



v i — v v, 



a 



- em) . 
en) . 



x i — Y em + (1 H u) [x 

/■> a \ i 
y i->- en + (1 u) [y - 

The metric induced on T m ^ n by this mapping is, 

ds 2 = 2dudv + ((1 + f u)dx + f (x - em)du) 2 + ((1 - ^u)dy - f (y - en)du)' 



(10) 
en by 

(11) 



(12) 



Together, the patches T m n cover the entire u > side of the grid. The metric ([12J) may 
be viewed as a metric on the union of all T m ^ n , R >0 x R 3 , if m and n are interpreted as 
functions of x and y respectively. That isQ 



m{x) 
n(y) 



[x/e + l/2\, 
[y/e + l/2\. 



(13) 



4 Here [x\ denotes the floor function, which returns the largest integer smaller than or equal to x. 
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These functions are not continuous, so we need to worry whether (12) describes a well- 
defined geometry on the whole u > half-space. The appropriate condition [131 E3] is t na t 
the metric induced on the common boundary of two patches is the same in both patches. 
On the x — e(m + |) boundary the metric on T m n induces the 3-dimensional metric, 

ds 2 = 2dudv + ((1 - f u)dy - f (y - en)du) 2 . (14) 

On the other side of the boundary the metric on the T m+ i n , can easily be seen to induce 
the same 3-dimensional metric. Similarly, we find that the 3-dimensional metrics induced 



on the y = e(n+ |) boundaries match, and hence that the metric (12) is well-defined on the 
whole u > half-space. 

Moreover, the 3-dimensional metric induced on the u = hyperplane is 

ds 2 = da; 2 + dy 2 , (15) 

which agrees with the metric induced from the Minkowski side. Consequently, the metric 
(|9| for u < and the metric ( 12 ) for u > together describe a well-defined geometry on the 



whole spacetime, given by the metric 

ds 2 = 2dudv + ((1 + fu6(u))dx + f (x - tm)9(u)du) 2 

+ ((1 - 1u6(u))dy -%y - en)9{u)du) 2 , (16) 

where 8(u) is the Heaviside stepfunction. 

This metric should describe the geometry of our piecewise flat configuration of defects. 
As a check, we can calculate the Riemann tensor of this metric. Since the metric is dis- 
continuous along the boundaries of the patches, we need to treat the curvature tensors as 
tensor distributions as described in [TJ]. The upshot of this is that in a case like this we 
can and should treat the metric components as distributions rather than normal functions. 
Using the fact that in a distributional sense 

A = S(x -i) = A(x), and (17) 

d 2 

u6{u) = 5{u), (18) 



du 2 

we find that the Riemann tensor is given by0 

R = 4a5(u)(A € (y + §)(du A dy) 2 - A e (x + §)(du A dx) 2 ), (19) 
where A e (x) is the Dirac comb with period e, i.e. 

A e {x) = ^6{x-ei). (20) 



We observe that the Riemann tensor indeed vanishes everywhere, except on the grid of 
defects where it has singular delta peaks as expected. We have therefore found a metric for 
our piecewise flat spacetime consisting of a grid of lightlike defects. 



5 For notational convenience we denote the tensor product of a 2-form with itself by a square, i.e. (du A 
dy) 2 = (du A dy) ® (du A dy). 
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Our next step is to find the continuum limit of this spacetime. For this we will send 
the grid size e to zero while keeping the density of defects (approximately) constant. The 
number of defects in a unit area of the grid scales inversely proportionately to e, hence the 
amplitude of the defects a should scale linearly with e to keep the total amplitude per unit 
area constant. Inserting a = ea^ in the above metric and taking the limit of e goes to zero, 
we find 

ds 2 = 2dudv + (1 + a u6(u)) 2 dx 2 + (1 - a u9(u)) 2 dy 2 , (21) 
where 9(u) is the Heaviside step function. The Riemann tensor (distribution) is given by 

R = Aa 5(u) ((du A dy) 2 - (du A da;) 2 ) , (22) 

from which the Einstein tensor may be obtained, which is seen to vanish. Consequently, 



the metric (21) is an exact solution of the vacuum Einstein equation, a gravitational wave 
consisting of a single plane wavefront. Waves of this type are called impulsive gravitational 
plane waves, and were first considered by Penrose [13] . 



IV. GENERAL GRAVITATIONAL PLANE WAVES 

We have succeeded in explicitly constructing a family of exact defect configurations that 
converges to a specific exact gravitational wave solution of general relativity in the limit that 
the configuration becomes continuous. In the construction we made some explicit choices 
for the sake of simplicity. We chose to align the grid with the coordinate axes and as a result 
described a plane wavefront with "+" polarization. Rotating the grid simply results in a 
rotation of the dx and dy terms in ( [2l| , producing a wavefront with a different polarization. 
For example, rotating the grid by 7r/4 produces a wavefront with "x" polarization. The 
result always is an impulsive plane wave with constant polarization. 

We could also take a different geometry for the grid. Equation ^ shows that the com- 
ponents of the grid do not need to be perpendicular, nor do the components have to have 
matching amplitudes. As long as the curvature density of the positive curvature defects is 
equal to the curvature density of the negative curvature defects, the Einstein tensor will 
vanish in the continuum limit. Different geometries will not produce any wavefronts that 
could not be produced by a square grid. 

Gravitational waves with longer wave packets may be constructed as a succession of 
wavefronts. Each wavefront contributes to the Riemann curvature in the following way, 

R = a5(u — Mo) [(dit A (cos0 dx + sin0 dy)) 2 — (du A (sin0 dx — cos0 dy)) 2 ] , (23) 

where <fi is the direction of the polarization of the wavefront. 

A continuous succession of wavefronts with a and 4> varying with u therefore produces a 
wave with curvature, 

R = a(u) [(du A (cos <p(u) dx + sin0(w) dy)) 2 — (du A (sm<p(u) dx — cos<p(u) dy)) 2 ] . (24) 

This is the most general curvature produced by exact gravitational plane waves of the type 
described by Bondi et al. in [T2] . 
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Figure 4. To first order, a piecewise flat cylindrical wave front should consist of axial (black) and 
azimuthal (white) defects moving inwards in the radial direction. The junction conditions imply 
that any additional radial defect (gray) should be present at each junction. 

V. OTHER GRAVITATIONAL WAVES 

We have shown that it is possible to construct general gravitational plane waves in our 
piecewise flat model of gravity. Since a general gravitational wave will locally resemble a 
plane wave, this suggests that it should be possible to construct any gravitational wave. This 
is further supported by the observations that any gravitational wave may be obtained in the 
linear limit of the model pQ|, and that its is possible to find a piecewise flat approximation to 
any solution of general relativity [IB]. Although this is all very suggestive, it is not enough 
to conclude that any gravitational wave may be constructed in our piecewise flat model of 
gravity. It is therefore instructive to explicitly examine the construction of a non-planar 
wave in our model and see what new issues arise. 

As an example, we will consider cylindrical gravitational waves. In the case of cylindrical 
symmetry the solutions of the vacuum Einstein equation can found by solving a single 
linear differential equationJTTl [T7] . Consequently, the complete set of possible solutions is 
known. In particular, it is known that it is impossible to form isolated impulsive wavefronts. 
Instead, impulsive wavefronts leave a wake of backscattered gravitational waves. Explicit 
solutions with impulsive wavefronts [TS] show that the form of this backscattered component 
can be quite complicated. This suggests that it may be too much to hope for a simple 
description of the backscattered component in terms of a piecewise flat geometry]^] Instead 
we will focus on constructing an impulsive cylindrical gravitational wave front, while we 
allow the backscattered component to be any defect configuration, i.e. not necessarily one 
with vanishing energy-momentum. 



6 This says nothing about whether a piecewise flat description exists. We simply observe that constructing 
one will be very complicated. 
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Let us consider a single ingoing cylindrical impulsive wavefront located at u = which 
propagates into a Minkowski background (i.e. the Riemann tensor vanishes for u < 0). In 
the limit of large radii (v — > oo), the wavefront will resemble an impulsive plane wave. Our 
first order guess for the piecewise flat cylindrical wavefront therefore consists of a grid of 
perpendicular massless defects with positive and negative energy oriented in the axial and 
azimuthal directions and moving in the radial direction (see figure [4]). 

Curving a flat grid to a cylinder introduces two new complications. The first is that due 
to simple geometric considerations the velocity of the azimuthal defects must be smaller than 
the velocity of the axial defects. The azimuthal defects, therefore, are not quite massless. 
However, as the grid is made finer this effect disappears and the azimuthal defects again 
become massless in the continuum limit. The second effect is that the four defects of the 
grid meeting at each junction do not lie in one plane. A consequence of this is that the 
junction condition ^ cannot be satisfied by these four defects. To satisfy (|5| a fifth defect 
pointing in the outward radial direction must be added to each junction (depicted as the 
gray lines in figure [5]). Just like general relativity, the piecewise flat model does not allow 
isolated cylindrical wavefronts. 



Repeating the steps of section III we can find the piecewise flat metric associated to this 



configuration and calculate the continuum limit. The result is 

,9 , , (u + v + au9(u)) 2 , , 9 . „. , . . 

ds 2 = 2dudv + ^ K -^d<f) 2 + (1 - au6(u)) 2 dz 2 . (25) 

The corresponding energy-momentum tensor is 

v — 1 1 

T = a 5(u)du 2 — 2a- — ; ——9(u)dudv. (26) 

v v ' (1 -au)(v+ (1 + a)u) y 1 v 1 

As expected the energy-momentum vanishes in the u < region, but does not in the 
'backscattered' u > region due to the presence of the radial defects. However, the energy- 
momentum also has a non-vanishing delta peak at u — for v ^ 1, i.e. it fails to model a 
cylindrical impulsive gravitational wavefront. This failure is the result of the two singular 
components of the Riemann tensor, 

R = Aa5(u) {du Adz) 2 - A-5{u) (du A -^=d0) 2 + 6{u) (•••), (27) 

v v 2 

scaling differently with v. The expected behaviour — as can be deduced from linear cylin- 
drical waves, or by expanding the exact cylindrical solution in [18] for small values of u - 
is for both curvature components to scale inversely with y/v. 

On the level of the piecewise flat geometry, this happens because the number of axial 
defects stays constant while the radius changes, causing their density to scale inversely with 
the radius. The density of the azimuthal defects, on the other hand, stays constant. To 
change the scaling of the impulsive front, the energy densities of the defects must change as 
the defects move inward. The only way this can happen is if the defects periodically emit 
new defects into the backscattered region. 

We implement this by having the ingoing impulsive front emit an outward moving lightlike 
front of perpendicular defects. A single cell of the piecewise flat geometry then takes shape 



7 We use cylindrical lightcone coordinates defined by u = and v = ^J-. 
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Figure 5. The geometry of a single cell of the piecewise flat geometry approximating an impulsive 
cylindrical gravitational wave front. 



shown in figure [5j By adjusting the frequency and the energies of the emitted perpendicular 
components we can achieve any scaling with v that we want. To achieve the proper scaling 
for the curvature of a cylindrical wave front, we choose the emitted components such that 
the metric in the continuum limit becomes, 



ds 2 = 2dudv + 



(u + v + cty/vu9(u)y 



u 



+ (l-a-^6(u)ydz 2 . 



The singular part of the corresponding energy-momentum tensor is 

r = o + 0(u)(--.). 



(28) 



(29) 



WThe impulsive wavefront has thus vanished from the energy-momentum tensor. Yet, the 
singular part of the Riemann curvature, 



R 



a 



[^=5(u)(du A dz) 2 - A^=5(u)(du A ^f<P? + 9{u) (••■), 



v 



(30) 



does not vanish, and has the expected scaling for a cylindrical impulsive wavefront. The 
described configuration indeed contains a cylindrical impulsive wavefront. 

The described piecewise flat geometry is just one way an impulsive cylindrical gravita- 
tional wavefront could evolve in our piecewise flat model of gravity. Recall that the model 
fundamentally describes gravity coupled to matter. As discussed in our previous articles 
[2j [3], the dynamics of the model is incomplete]^] It needs to be complemented by addi- 
tional rules which relate hoe the matter content interacts, much like general relativity needs 
to be complemented by equations of motion for the matter content. Different dynamical 
completions of the model lead to different evolutions of cylindrical wavefronts. The specific 



This is related to the different ways in which a piecewise flat geometry can be continued after the collision 
of two defects. 



12 



evolution described here corresponds a matter interaction where interacting gravitational 
waves can produce other types of matter. 

One might ask if there exists a dynamical completion of the model that corresponds to 
pure gravity with no matter present at all. Of course, this is of academic interest only since 
the physical world also contains matter. In such a completion the backscattered compo- 
nent of the evolving cylindrical wavefront must also satisfy the vacuum Einstein equation. 
This requires the addition of additional defects to the configuration. These will break the 
remaining residual symmetry of the piecewise flat geometry. As a consequence it will no 
longer be possible to describe the geometry as an infinite collection of identical piecewise flat 
cells, substantially adding to the complexity of the description. This breaking of symmetry 
indicates that cylindrical gravitational waves may not be the best setting to further study 
the dynamical completion of the piecewise flat model. 

Nonetheless, the study of the interaction of gravitational waves in the piecewise flat 
model may provide essential guiding insights for dynamically completing the model. A 
more promising setting for this study is the interaction of two plane wavefronts. This 
problem is well documented in the case of general relativity coupled to various types of 
matter [T9l - |2~2"] . Moreover, the intrinsic plane symmetry is much more suited to a piecewise 
flat geometry. Furthermore, it is well known that the future of plane wave collisions can 
contain regions which are locally diffeomorphic to black hole solutions [121 [22]- Studying 
plane wave collisions may therefore offer crucial insight in one of the major open problems 
of the piecewise flat model; the construction of non-trivial stationary vacuum solutions such 
as black holes. However, the study of dynamical completions of the model is a major project 
on its own and will not be attempted here. 

VI. CONCLUSIONS 

In this article we have extended our previous linear result p] for obtaining gravitational 
waves in our piecewise flat model for gravity in 3+1 dimensions to an explicit construction 
of a family of exact piecewise flat configurations that approaches exact plane wave solutions 
of general relativity. This suggests that the construction of general gravitational waves may 
be possible. As an example, we have constructed the evolution of a cylindrical gravitational 
wavefront allowing the backscattering to produce other types of matter. 

Understanding the construction of gravitational waves is an important step towards show- 
ing that the proposed model can indeed serve as description of gravity in 3+1 dimensions. 
Further study of the interactions of gravitational waves may provide essential guiding in- 
sights for completing the incomplete dynamics of the model. 
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